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The double logarithmic terms a s In 2 x are important to predict precisely the small 
x behavior of the spin structure function gi . We numerically analyze the evolution 
of the flavor non-singlet gi including the all-order resummed effect of these terms. 
It is pointed out that the next-to-leading logarithmic corrections produce an un- 
expectedly large suppression factor over the experimentally accessible range of x 
and Q 2 . This implies that the next-to-leading logarithmic contributions are very 
important in order to obtain a definite prediction. 



1 Introduction 

Recently many experimental and theoretical works have been devoted to the 
polarized structure function g\. Especially, it is very important and desirable 
to know the small x behavior of 171 in the light of the Bjorken and Ellis- Jaffc 
sum rules. Since the verification of these sum rules requires the knowledge of 
the structure function over the entire x regions, one has to rely on the theoret- 
ical prediction in the experimentally inaccessible small x region. Although the 
Regge prediction has usually been assumed for the extrapolation of the experi- 
mental data to the small x region, the recent data show a clear departure from 
this naive Regge prediction El. This fact means that the perturbative QCD 
effects are very important. 

Recently, various extrapolations have been proposed using the DGLAP equa- 
tion El. However, it has been known that there appear the double logarithmic 
terms a s ln 2 x in the perturbative calculations. These logarithmic corrections 
seem to give large effects in the small x region and are important to get more 
reliable-predictions on the small x behavior of g%. Bartels, Ermolaev and 
Ryskintj have given the resummed expression for the partonic structure func- 
tion gP arton D y using the Infra-Red Evolution Equation and confirmed the old 
result by Kirschner and Lipatovcl. They claim that the resummed effects may 
be important. But, when extracting the physical structure function of hadrons 
from the partonic one, there is possibility that their conclusion at the parton 
level is not necessarily true. Indeed, the recent numerical analysis by Bliimlein 
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and Vogtu shows that there are no significant contributions from the resum- 
mation of the leading logarithmic (LL) corrections in the HERA kinematical 
region. The different conclusions between at the partonic and hadronic level 
might come fronLthe fact that the resummed "coefficient function" was not 
included in Ref. □ because it falls in the next-to-leading logarithmic (NLL) 
corrections and depends on the factorization scheme adopted. It is also to be 
noted that the slightly steep input density was used in Ref. til. The evolution, 
in general, strongly depends on the input parton densities. If one chooses a 
steep input function, the perturbative contribution will be completely washed 
away. So it will be also interesting to see the sensitivity of the results to the 
choice of the input densities. 

In this report, we reanalyze the numerical impact of the resummed effects on 
the small x behavior of g\ (non-singlet part). The coefficient function can not 
be included consistently at present since the anomalous dimension has been 
calculated only to the LL order. However we consider the effects of this part 
because we could firstly clarify an origin of the above different conclusions and 
secondly get some idea about the magnitude of the NLL order corrections in 
the resummation approach. We also consider two different input densities: one 
is flat corresponding to the naive Regge prediction and the other is. steep in 
the small x region. Details of the calculations may be found in Ref. B. 



2 Resummed structure function gi 

The flavor non-singlet part of g\ in the moment space is given by, 

91 (Q 2 ,N) ee J dxx N - 1 g 1 {Q\x) = ^C(a s (Q 2 ), N)Aq(Q 2 , N) , 
where 

C(a s (Q 2 ),N) ee f dxx N - l C{a s {Q 2 ),x) , Aq(Q 2 7 N) = f dxx N - x Aq(Q 2 , x) , 
Jo Jo 

and Aq(Q 2 , x) (C) is the flavor non-singlet combination of the polarized parton 
densities (the coefficient function), (e 2 ) is the average of the quark's electric 
charge. The DGLAP equation is, 

Q 2 jL-Aq{Q 2 , N) = - 1 {a s {Q 2 ),N)Aq(Q 2 , N) . (1) 

Here the anomalous dimension 7 is the moment of the "splitting" function. 
The coefficient function C(a s ,N) and the anomalous dimension j(a s ,N) can 
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be calculated perturbatively and are expanded in the powers of a s , 



C(a s ,N) = l + J2c k (N)a': , ~/(a s , N) =^T 7 k (N)a k s . 

k=l k=l 

where a s = . The singular behaviors of the anomalous dimension and 
the coefficient function as x — > appear as the pole singularities at N = 0. 
We must resum these singular terms to alLorders to get a reliable prediction. 
This resummation has been done in Refs. Uu and the resummed part of these 
functions are given as follows, 

j(a s ,N) = Urn >y(a B ,N) = -fo{N)/^ 2 , (2) 

N 

= lim C(a s ,N) = r- . (3) 

N - f^ } (N)/8ir 2 




a s 



2tt 2 N J8 



16ir 2 N c a s 4cMe z2/A D-i/2N*(z)) 



D p {z) is the parabolic cylinder function. One can easily see by expanding 
Eqs.(||j3|) in terms of a s that the resummed expressions Eqs.(||,||) reproduce 
the known NLO results in the MS scheme. Therefore, it is quite plausible that 
Eqs.(||j3|) correctly sum up the "leading" singularities to all orders. 

Before going to the numerical analysis, we explain the reason why the re- 
summed coefficient part was discarded in the analysis of Blumelein and Vogt. 
For definiteness, let us use the so-called DIS scheme. The parton density and 
the anomalous dimension in the DIS scheme are obtained by making the trans- 
formations, 

Aq -> Aq DIS = CAq , 7 DIS = C jC- 1 - /?(a.) In C. (4) 

oa s 

Using the resummed 7 and C Eqs.(|],[|), we get the anomalous dimension in 
the DIS scheme, 

k / / ri, \ k^ 



^-VIM^) + + J*(*) , (5) 



k=l k=2 
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where the second term comes from the resummed coefficient function and d 
are numerical numbers independent of N . Since the second term has an extra 
N in comparison with the first term in the each power of a s , the resummed 
coefficient function belongs to the NLL order corrections in the context of the 
resummation approach. This is the reason why the LL analysis of Blumlein 
and Vogt dropped the resummed coefficient function. One must include the 
NLL order anomalous dimension, which has not yet been available, to perform 
a consistent analysis at the NLL level. 

3 Numerical consideration 

We show our numerical results. From Eq.(|j), in DIS scheme, g\ is given by 
Si M 2 ) = ^ ^- N ^ DIS (N,Q 2 ) (6) 

2 Jc-ioo 27T l 

The DGLAP equation Eq.([j]) is easily solved with anomalous dimension j DIS , 
Aq DIS (N,Q 2 ) = exp (- f*^ A ^ (iV ,Qg). 

V J as (Ql) P J 

The anomalous dimension j DIS which includes the resummed effects is ob- 
tained by using Eqs.(§,| ) and the relation Eq.(|), 

7 ms (iV) = a sl \N)+a 2 sl 2 (N)+K(N, In (l + a s c l + H(N, a s )) , 

da s 

(7) 

where 7 ' and c 1 are respectively the exact anomalous dimension and co- 
efficient function at the one and two-loop fixed order perturbation theory. 
K(N,a s ) (H(N, a s )) is the resummed anomalous dimension Eq.(^) (the co- 
efficient function Eq.(^)) with k = 1,2 (k — 0,1) terms being subtracted to 
avoid the double counting. Although the anomalous dimension at JV = 1 
should vanish due to the conservation of the (non-singlet) axial vector cur- 
rent, the resummation of the leading singularities in N doesjiot respect this 
symmetry. One of prescriptions [] to restore this symmetry isH, 

K(N,a s ) -+K(N,a.)(l-N) . 

In order to estimate gi , we must assume the appropriate function, for the input 
density Aq DIS (N, Qq). The explicit parameterization we use isl3, 

Aq DIS (Q 2 ,x)=N(a,P,a)r 1 x a (l-xf(l + ax) , 
''Our final conclusion remains the same qualitatively if we choose other prescription. 
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where N is a normalization factor such that J dxNx a (l — x)@(l + ax) = 1 and 
V = §9A (§A = 1-26) in accordance with the Bjorken sum rule and we choose 
the input scale Qq — 4GeV 2 . Note that the small x behavior is controlled by 
the parameter a. Since we are also interested in the sensitivity of the final 
results to the small x behavior of input densities, we choose two types of the 
input densities A and B: A is a function which is flat at small x (x a , a ~ 0), 
and B is slightly steep (a ~ —0.2) . A and B correspond to the following values 
of parameters, 

A (B) : a = +0.0 (-0.2) , (3 = 3.09 (3.15) , a = 2.23 (2.72) . 

We put the flavor number n/ = 4 and Aqcd = 0.23GeV. 

Now, let us explain how to perform the Mellin inversion Eq. (j6j) which is an 
integral in the complex TV-plane. The contour integration along the imagi- 
nary axis from c — ioo to c + ioo is numerically inconvenient due to the slow 
convergence of the integral in the large \N\ region. To get rid of this prob- 
lem, we deformed the contour to the line which have an angle <fi (cf> > n/2) 
from the real N axis. By this change of the contour, we have a damping fac- 
tor exp(|7V| ln(l/a;) cos</)) which strongly suppresses the contribution from the 
large \N\ region. In the integration along this new contour, we will be able to 
cut the large \N\ region. Finally we have checked the stability of results by 
changing the contour parameter. One can find the details of this technique in 
Ref.H. 

First we estimate the case which includes only the LL correction. The evo- 
lution kernel in this case is obtained by dropping H(N,a s ) in Eq.(|7|). Fig. la 
(lb) shows the LL results (dashed curves) after evolving to Q 2 = 10, 10 2 , 
l0 4 GeV 2 from the A (B) input density (dot-dashed line). The solid curves are 
the predictions of the NLO-DGLAP evolution. These results show a tiny en- 
hancement compared with the NLO-DGLAP analysis and are consistent with 
those in Ref.El The enhancement is, as expected, bigger when the input density 
is natter. However any significant differences are not seen between the results 
from different input densities. Next, we include the NLL corrections from the 
resummed "coefficient function". We show the results in Fig. 2 by the dashed 
curves. The results are rather surprising. The inclusion of the coefficient func- 
tion leads to a strong suppression on the evolution of the structure function at 
small x. To understand our numerical results, it will be helpful to remember 
the perturbative expansion of the resummed results Eqs.(H||). Using the ex- 
plicit values Nc = 3, Cf = 4/3, we obtain for the anomalous dimension in the 
DIS scheme Eq.(|), 



jDIS = N 
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Figure 1: The LL evolution as compared to the DGLAP results with the flat input A (la) 
and steep one B (lb). 
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Here note that: (1) the perturbative coefficients of the LL terms (the first 
part of Eq.(ra)) are negative and those of the higher orders are rather small 
number. This implies that the LL corrections push up the structure function 
compared to the fixed-order DGLAP evolution, but the deviations are expected 
to be small. (2) the perturbative ones from the NLL terms (the second part 
of Eq.(||)), however, are positive and somehow large compared with those of 
the LL terms. This positivity of the NLL terms has the effect of decreasing 
the structure function. It might be also helpful to assume that the saddle- 
point dominates the Mellin inversion Eq.(^). We have numerically estimated 
the approximate position of the saddle-point and found that the saddle-point 
stays around N$p ~ 0.31 in the region of x ~ to 10 -2 . By looking 
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Fig.2a Fig.2b 

Figure 2: The NLL evolution as compared to the DGLAP results with the flat input A (2a) 
and steep one B (2b). 

at the explicit values of the coefficients in Eq.(||), the position of the saddle- 
point seems to suggest that the NLL terms can not be neglected. Since the 
coefficients from the higher order terms are not so large numerically, it is also 
expected that the terms which lead to sizable effects on the evolution may 
be only first few terms in the perturbative series in the region of x we are 
interested in. Fig. 3a (3b) shows the numerical results of the contribution from 
each terms of the NLL corrections in Eq.((|) at Q 2 = 10 2 GeV 2 with the A (B) 
type input density. The solid (dot-dashed) line corresponds to the NLL (LL) 
result. The long-dashed, dashed and dotted lines correspond respectively to 
the case in which the terms up to the order a 2 , a^, af, are kept in the NLL 
contributions. One can see that the dotted line already coincides with the full 
NLL (solid) line. These considerations could help us to understand why the 
NLL corrections turns out to give large effects. 

4 summary 

We have numerically studied the small x behavior of the flavor non-singlet gi 
by taking into account the resummed effect of a s ln 2 x. Our LL analysis is con- 
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Figure 3: Contributions from the fixed order terms in the NLL resummation with the flat 
input A (3a) and steep one B (3b). 



sistent with the results by Blumlein and VogtH. We have also performed the 
analysis which includes a part of the NLL corrections from the resummed coef- 
ficient function in the light of the assertion of Bartels, Ermolaev and Ryskind, 
though this is not theoretically consistent. Our results suggest that the LL 
analysis is unstable in the sense that a large suppression effect comes from the 
resummed coefficient function which should be NLL corrections. We need a 
full NLL analysis to make a definite conclusion. 

This work was supported in part by the Monbusho Grant-in- Aid Scientific 
Research No. A (1) 08304024, No. C (1) 09640364. 
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